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We study transport in ferromagnetic-superconductor/normal- metal systems. It is shown that 
charge and spin currents are pumped from ferromagnetic superconductors into adjacent normal 
metals by adiabatic changes in the order parameters induced by external electromagnetic fields. Spin 
and charge pumping identify the symmetry of the superconducting order parameter, e.g., singlet 
pairing or triplet pairing with opposite or equal spin pairing. Consequences for ferromagnetic- 
resonance experiments are discussed. 
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Ferromagnetism induces a spin-dependent asymmetry 
in the densities of itinerant carriers. In contrast, su- 
perconductivity pairs electrons with equal or opposite 
spins depending on the symmetry of its order parameter. 
The coexistence of the two order parameters has been 
considered to be a rare phenomenon. However, recent 
experimental progress has demonstrated that ferromag- 
netism and superconductivity coexist in some materials 
like RuSrsGdCusOg Q, UGes 0, ZrZnj 0, and URhGe 
The experiments find triplet pairing in URhGe, and 
strong indications of triplet pairing in UGe2 and ZrZn2; 
they furthermore suggest that the same electrons are re- 
sponsible for ferromagnetism as well as superconductiv- 
ity. Besides, ferromagnetism and superconductivity are 
predicted to be simultaneously induced in hybrid ferro- 
magnet (F)/normal-metal (N)/superconductor (S) sys- 
tems j5|- These recent experimental discoveries, and 
the possibility of tailoring superconductivity and ferro- 
magnetism in nanoscale systems, enable exploring novel 
physics involving pairing and spin-related transport pro- 
cesses. A variety of interesting spin phenomena have al- 
ready been observed in hybrid F/N and semiconductor 
systems: e.g., giant- magnetoresistance — related effects, 
spin precession, and current-induced magnetization dy- 
namics . It is therefore natural to expect that interest- 
ing rich phenomena should also occur in FS/N systems. 

This Letter demonstrates how the coexistence of su- 
perconductivity and ferromagnetism is manifested in the 
adiabatic pumping in hybrid FS/N structures. In par- 
ticular, we study spin and charge pumping when the 
magnetization slowly precesses, which can be achieved in 
ferromagnetic resonance (FMR) and in current-induced 
magnetization dynamics. FMR experiments have al- 
ready been carried out to investigate the magnetism in 
RuSr2GdCu208 0- We also consider pumping by slow 
variations in the phase of the singlet or triplet order pa- 
rameter, or in the direction of the triplet order parameter. 
These can be induced by electric and magnetic fields to 
be discussed below. By pumping we thus mean the spin 
and charge flows into the adjacent normal conductors in 
response to adiabatic changes in the FS order parameters. 



Consequently, in the case of pumping by the magnetiza- 
tion direction, we compute the spin current Is and the 
charge current Ic for a given rate of the magnetization- 
direction change, and for pumping by changing pairing, 
we compute the same quantities as functions of the phase- 
change or the direction-change rates of the pair correla- 
tions. We employ two approaches giving identical re- 
sults: 1) solving the time-dependent ac problem directly 
and 2) using a gauge transformation to obtain a time- 
independent dc problem. We first explain our model and 
notation, before proceeding to the derivation and results. 
Experimental consequences are discussed in the end. 

A ferromagnetic superconductor is treated in the 
mean-field approximation, where ferromagnetism is rep- 
resented by the average magnetization and superconduc- 
tivity is described by a pair potential. Our model is phe- 
nomenological and we do not discuss the microscopic ori- 
gin of the exchange field or the superconducting pairing. 
The Bogoliubov-de Gennes (BdG) equation is 
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where ^ = H{)\ + a ■ Cxc is the single-particle Hamilto- 
nian, e^c = Cxcin (with exc > 0) is the ferromagnetic 
exchange field along the magnetization direction m and 
^ — {dol + a ■ ctj i&y is the superconducting pair po- 
tential, in terms of the singlet (scalar) part, do, and the 
triplet (vector) part, d — {dx,dy,dz)- = {v) ju'^) are 
spin-dependent electron wave functions and iP" = {v\ v-^) 
are those of holes; = {cFx, Oy^a^') is the vector of Pauli 
matrices. The single-particle Hamiltonian Hq contains 
the kinetic and potential-energy terms. The Fermi en- 
ergy is taken to be the largest relevant energy scale. The 
local exchange field can be position dependent close to 
the interface, exc(r), and the pair potential A(k, r) can 
also be position as well as wave- vector, k, dependent 0. 
For simplicity, we assume that the exchange field and the 
pair potential are uniform inside the superconductor and 
drop to zero at the FS/N interface: exc(i") — exc6(z) 
and A(k, r) — A(k)0(z), where 8(2:) is the Heaviside 
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step function and z is the coordinate perpendicular to 
the FS/N interface. 

Fermionic statistics dictates A(k) — — A-'"(— k) 0. 
The singlet (triplet) part of A thus needs to have even 
(odd) parity: do(k) = do(-k) and d(k) = -d(-k). We 
study in the following two simple cases of triplet super- 
conductors: opposite-spin pairing (OSP) along the ex- 
change field, d(k)x Cxc = and equal-spin pairing (ESP), 
d(k) • — [lOj- Triplet OSP superconductors are 
described by a (complex- valued) vector d(k) — (i(k)m 
along the magnetization direction. We show that the 
transport properties in triplet OSP are similar to those 
of singlet pairing. In ESP, superconductivity occurs inde- 
pendently for spins along and opposite to the magnetiza- 
tion direction: By choosing the magnetization along the z 
axis, the superconducting pair potential decomposes into 
two terms, corresponding to spins up and down along the 
z axis, S^^^i^ = =Frfx(k) -l-idy(k). Since superconduct- 
ing correlations do not mix the spin-up and down sub- 
systems, two ESP phases can be distinguished: the Af 
phase, where pairing occurs only for spin a (i.e., (P ^ Q 
and dr" = 0) and the A2 phase, where pairing occurs for 
both spins. A large exchange interaction, Cxc ^ If^ol (Ml): 
suppresses superconducting singlet (triplet) OSP corre- 
lations 0. We therefore only consider Exc < Mo I (Ml) 
for these systems, so that the quasiparticle excitations 
have a finite gap. Triplet ESP have a quasiparticle gap 
in the superconducting spin channels independent of the 
size of the exchange interaction [ij3 | , and we do not make 
assumptions about the ratio of the exchange field to the 
pair potential in this case. 

Let u s apply the standard scattering-matrix approach 
[HElElElEi to an FS/N system. We assume that 
the Hamiltonian is continuous across the FS/N in- 
terface and incorporate interfacial disorder and band- 
structure mismatch into a "disordered" normal region 
[llj |. Similarly to Ref. 0|, we solve the BdG equation 
for an electron (or hole) incident on the specular FS/N 
interface from the normal-metal side. The total reflec- 
tion matrix is then found by concatenating the FS/N 
reflection with the scattering by the normal disordered 
region for electrons and holes. The FS layer, in series 
with the disordered region, is viewed as a scatterer for 
electrons supplied by the normal reservoir, see Fig. ^ 
The problem is simplified in the clean-superconductor 
limit, where the FS mean free path is longer than the su- 
perconducting coherence length HvfKt^^), expressed in 
terms of the Fermi velocity vp and quasiparticle gap A. 
Focusing on the low-temperature regime, fc^T ^ A, we 
define the M x M spin-dependent electron— ^electron and 
electron^hole reflection matrices and rj^^, where M 
is the total number of quantum channels, or transverse 
waveguide modes, at the Fermi level of the normal-metal 
lead and a is the spin label along the magnetization di- 
rection for incident electrons. For singlet or triplet OSP, 
rf describes reflection into holes with spin — cr, while 
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FIG. 1: (a) A ferromagnetic superconductor coupled to a nor- 
mal reservoir through a specular FS/N interface in series with 
a normal disordered region, (b) Electrons incident on the 
FS/N interface from the normal side are reflected as electrons 
or holes and transmitted as decaying electron- or hole-like 
quasiparticles. 



for triplet ESP, reflected holes have the same spin a as 
incident electrons. 

We find that for both dc and adiabatic ac transport, 
spin and charge currents are uniquely determined by 
two real-valued conductances and one complex-valued 
conductance, representing transport of spins aligned to 
the magnetization, antialigned to the magnetization, and 
transverse to the magnetization, respectively: 
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It is convenient to define the total conductance g = + 
g^^ and the polarization p = {g^^ ~9^^)/ 9- We will below 
interpret how these conductances determine charge and 
spin flow, and discuss their values in various limits. 

Let us first fix the magnetization direction m and con- 
sider pumping by an adiabatically varying phase 0(t) 
of the superconducting pair potential, dt4> — uj 
(exc, A)/?i. As is well known, this time-dependent prob- 
lem can be transformed into a dc problem by a gauge 
transformation U{t) = e*"^*^*^/^: Varying the phase re- 
sults in nonequilibrium transport corresponding to a dc 
response at a voltage bias hioji^e). The pumped cur- 
rent is thus nothing more than a response to a voltage 
V = huj/{2e) apphed between the superconducting con- 
densate and the normal metal, and we compute 
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g is thus the usual Andreev-reflection conductance, in 
units of /{ttK). The spin current is determined by the 
conductance polarization p. Eqs. I@J can also be derived 
in the ac pumping framework of Ref. (I5| . 
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Pure spin flow is generated by the variations in the 
magnetization direction m, induced by, e.g., a resonant 
rf magnetic field. In Ref. we computed the corre- 
sponding pumped spin flow into the normal metal, 

h { ^, dm dm.\ , . 

Is = — Reg'^m X h Img'-^ , (5) 

47r \ dt ' dt J 

for a general F/N contact with no superconducting cor- 
relations. Here g"^^ = Tr[l — ?'|e(''ie)^] F/N mix- 
ing conductance for transverse spins expressed in terms 
of the normal-side reflection matrices |l3l |. The pumped 
spin flow induces an enhanced Gilbert damping when the 
normal metal is a good spin sink, so that g'^^ is exper- 
imentally measurable |3- We have generalized Eq. (O 
to spin pumping by a superconducting ferromagnet: The 
pumped spin flow remains of the form ((SJ with a redeflned 
mixing conductance g^^ (jSJ. There is no accompanying 
charge pumping. We derive Eqs. Q and ij^l by two differ- 
ent methods: First, we extend the pumping approach of 
Ref. 01 , where the scattering matrix is time-dependent 
in spin space due to a slow variation in the magnetization 
direction m, to include electron^hole reflection. Sec- 
ondly, and much simpler, the calculation is reformulated 
as a dc problem in the spin frame of reference which is 
moving together with m.{t). The latter is achieved by (in- 
stantaneously) applying a spin-rotation operator around 
the vector m x 9tm and correspondingly adding a new 
term in the Hamiltonian: ^' = — (/i/2)(m x 9tm) • a. 
This term corresponds to an equilibrium transverse spin 
accumulation which, in turn, induces the spin current 
(O. The mixing conductance Q is obtained after ex- 
tending the F/N dc theory of Ref. 13] to account for 
electron^hole reflection. By unitarity of the scattering 
matrix the real part of g^^ is bounded from above 
by twice the number of channels: Reg^-'- < 2M. 

For triplet OSP superconductors, the derivation lead- 
ing to Eq. (jHJ assumes that the triplet pair-potential di- 
rection moves together with m, while its wave-vector 
dependence is locked to the atomic lattice: d(k, t) = 
(i(k)m(t). For triplet ESP, the situation is more complex 
since the vector pair potential d is perpendicular to m, 
resulting in an additional dynamic degree of freedom. In 
particular, the magnetization motion does not uniquely 
determine the trajectory of d. In deriving Eq. ((SJ), we as- 
sumed that d rotates together with m around m's instan- 
taneous rotation axis. We can get more complex trajec- 
tories by combining (instantaneous) rotations of m and 
d with "twisting" of d around m and the overall phase 
variation of the pair potential. The induced currents will 
then be given by adding the corresponding contributions 
to pumping. Calculating the exact trajectories for real- 
istic systems, which might be governed by spin-orbit in- 
teractions in the lattice field or other microscopic details, 
lies beyond the scope of this paper. A simple example is 
the clockwise rotation of d around m with frequency lo 



which induces additional currents 

e h 
Ic = -—pguj and Is = —guncn. (6) 
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Note that while the preceding equations are general, 
Eq. applies to the case of triplet ESP only. For a 
FS/N system disconnected from an Ohmic circuitry, the 
low- frequency charge flow vanishes, so that the overall 
phase variation of d and its twisting around m must re- 
sult in canceling charge currents, /c's, in Eqs. Q and 
(jnj, but a flnite net spin current for p < 1. Eqs. © 
can be derived similarly to Eq. lO, either as a time- 
dependent pumping problem or a dc problem in the 
gauge-transformed frame of reference rotating with d. 

Eqs. (0J, |5j, and 10) are general expressions for 
pumped charge and spin flows by varying FS order 
parameters. These currents are all governed by two 
real-valued and one complex conductance parameters: 
g^\ g^^, and g^^, which can be evaluated in micro- 
scopic models. A finite value of the electron— >hole re- 
flection coefficient, r^^, requires that an electron inci- 
dent from the normal-metal reservoir gets transmitted 
through the disordered region, converted into holes at 
the interface, and transmitted back into the normal lead 
as a hole, see Fig. ^ For small (characteristic) trans- 
mission eigenvalues T of the normal disordered region, 
g^^ and g^^ therefore scale as rather than as T, 
as in the Landauer-Biittiker formula for an N/N con- 
tact. Indeed, the Andreev conductance for a singlet 
nonmagnetic-superconductor/normal- metal contact was 
shown to be gs/N = 2?'m/(2 - T^)'^, where m la- 
bels the transmission eigenvalues for scattering in the 
normal metal the limit of no disorder (and no 

band-structure mismatch), Tm = 1 and (?s/n = 2M, i.e., 
Andreev reflection causes a doubling of the conductance, 
as compared to the N/N interface. We generalize gs/N 
to the case of a magnetic singlet superconductor with an 
s-wave symmetry of the pair potential, c?o(k) = |c?o|e"^: 

g = V — (7) 

^ Z.(2_r„)2-4(l-T„)(6WMo|)2 ^' 

and p — 0. The mixing conductance Q of s-wave mag- 
netic superconductors can also be expressed in terms of 
the scattering matrix of the disordered region, by gen- 
eralizing the formalism of Ref. j^. For that we need 
to concatenate the transfer matrix of the normal disor- 
dered region with electron<-!-hole conversion at the inter- 
face: r\,, , = e*('3±0) J i ^ g-«(/3T0) where 

eh{he} ehihe) ' 

(3 — arccos(exc/Mo|)- The situation is even more compli- 
cated for the k-dependent pair potential (which is always 
the case for the odd-parity triplet pairing) and we do not 
pursue it here. The exception is the case with no dis- 
ordered region in our model. Fig. ^ The conductance 
parameters for the singlet and triplet OSP in such bal- 
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listic systems are then given by 
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where the ± sign corresponds to singlet (triplet) pair- 
ing and A = |do| or \d\, respectively, k is the incident- 
electron wave vector at the Fermi level, |k| = kp, and 
kj_ is its transversal projection in the lead with cross 
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J(f}-<l>) and rL 



(3 = arccos(exc/M|) for the triplet OSP with d = Idle*"*, 
for a given k. 

It is instructive to discuss the mixing conductance g^^ 
(O values in some special cases of a clean interface with 
matched band structures since it can be measured exper- 
imentally in FMR. For an F/N interface, r^^ — r^^ — 0, 
while g^^ = M is large determined by the number of 
transverse wave-guide channels M, assuming that the F 
layer is thicker than the ferromagnetic coherence length 
hvp/{TTexc)- For a perfect electron^hole reflection off the 
singlet superconductor, we find 
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the limit when e^c A, resulting in a vanishing mix- 
ing conductance, as follows from Eq. (jS)). This is easy to 
understand since rj^^ = —rj^^ means that the transverse 
spin-t electrons get reflected as the spin-| holes which ex- 
actly cancel the incident transverse spin current. In the 
analogous limit of the triplet OSP, 
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doubling the F/N mixing conductance. We thus find that 
gTi = M, 0, 2M for F/N, singlet FS/N, and triplet OSP 
FS/N interfaces. Since g^^ is a direct measure of the fer- 
romagnetic Gilbert-damping enhancement [l^ , the onset 
and the nature of the superconducting pairing has non- 
trivial consequences in FMR experiments 7] . It is worth- 
while noting that g^-^ 2M is a rather special limit for 
the mixing conductance in multilayer magnetoelectronic 
circuit theories 0, For example, in the case of a 
symmetric FS/N/FS structure, it should result in a dy- 
namic locking of the two magnetizations which can be 
measured experimentally [l7| |. 

Finally, in triplet ESP with no disordered region. 
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in the A^ phase, and 
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= A 



{1 + exp (z[0T(k)-<^i(k)])}(12) 



in the A2 phase, 
ing spin channel with d"^ = \d° 
deriving Eq. Q , we made a convention for the instanta- 
neous coordinate system of x oc m x dtm, y cx —dtm, 



for each superconduct- 
^ 0. Note that in 



and z = m. [It is necessary to specify the coordinate- 
system convention in the case of the A2 phase, because 
i^T(i) _ zpdx + idy and, since d transforms as a vector, 
the relative phase of d^^-'-^ entering Eq. (|12|l depends on 
the choice of the x and y axes.] It then follows that the 
second term in the curly brackets of Eq. (|12|l is modu- 
lated for a small-angle precession of m, unless d twists 
around m with its precession frequency (i.e., following 
the instantaneous rotation axis). In the former case, the 
theory therefore predicts an anisotropic Gilbert-damping 
parameter. For the triplet ESP A2 phase, the symme- 
try of the superconducting order parameter may thus be 
seen in an additional anisotropic FMR line width, as one 
changes the temperature across the F-to-FS transition. 

In summary, we have studied the interplay between fer- 
romagnetism and superconductivity in adiabatic pump- 
ing by varying order parameters in FS materials in con- 
tact with normal metals. We demonstrate that the sym- 
metry of the superconducting pair potential is reflected 
in the conductance parameters which govern the pumped 
spin and charge flows. An experimental quantity of a par- 
ticular interest, which encodes information about both 
ferromagnetic and superconducting correlations, is the 
mixing conductance g^^ which governs the Gilbert damp- 
ing of the magnetization dynamics Consequently, 
signatures of the FS order parameter can be measured in 
thin film FS/N or FS/N/FS FMR experiments. 
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